Magnetocrystalline anisotropy and uniaxiality of MnAs/GaAs(100) films 



J. Magnus Wikberg-'^, Ronny Knut^, Sumanta Bhandary^, Igor di Marco^, Mikacl Ottosson"^, Janusz Sadowski"', 
Biplab SanyaP, Pal Palmgren^, Chcuk W. Tai'^, 011c Eriksson^, Olof Karis^ and Peter Svedlindh^ 
^Department of Engineering Sciences, Uppsala University, P.O. Box 534, SE-751 21 Uppsala, Sweden 



P.O. 

^Department of Physics and Astronomy, Uppsala University, Box 516, 751 20 Uppsala, 



Sweden 
Sweden 



Department of Materials Chemistry, Uppsala University, P.O. Box 538, SE-751 21 Uppsala. 
"^MAX-lab, Lund University, P.O. Box 118, SE-221 00 Lund, Sweden and 

'^Department of Materials and Environmental Chemistry, 
Arrhenius Laboratory, Stockholm University, SE- 106 91 Stockholm, Sweden 

We present an investigation of the magnetic behavior of epitaxial MnAs films grown on GaAs(lOO). 
We address the dependence of the magnetic moment, ferromagnetic transition temperature (Tc) 
and magnetocrystalline anisotropy constants on epitaxial conditions. From thorough structural 
and magnetic investigations, our findings indicate a more complex relationship between strain and 
magnetic properties in MnAs films than a simple stretch/compression of the unit cell axes. While 
a small increase is seen in the anisotropy constants the enhancement of the magnetic moment at 
saturation is significant. X-ray magnetic circular dichroism results show a behavior of the spin- 
and orbital-moment which is consistent with a structural transition at Tc. In particular, we find 
that the ratio of the orbital to spin moment shows a marked increase in the coexistence region of 
the ferromagnetic a- and paramagnetic /3-phases, a result that is well in accord with the observed 
increase of the c/a-ratio in the same temperature region. The ah initio density functional calculations 
reveal that the magnetic properties are more sensitive towards change in in-plane axis as compared 
to a change of the out-of-plane axis, which is explained by the analysis of band structures. The 
effects of electron correlation in MnAs using ab initio dynamical mean field theory are also presented. 

PACS numbers: 75.30.Gw, 75.70.-i, 71. 15. Mb 



I. INTRODUCTION 



The search for suitable materials for spintronic appli- 
cation has spurred a large interest in materials that can 
be deposited on ordinary semiconductors such as Si and 
GaAs. Bulk MnAs exhibits a ferromagnetic (FM) tran- 
sition temperature (Tc) of 313 K»i where it is accom- 
panied by a structural phase transition going from the 
FM a-phase to the orthorhombic /3-phase and as the 
temperature is raised to 398 K MnAs exhibits an ad- 
ditional phase transition back to a hexagonal 7-phase. 
The long term interest in MnAs lies in the first order 
magnetostructural phase transition with a Tc that can 
be altered by for instance applying a stress to the mate- 
rial resulting in tensile/comprcssive strain of the MnAs 
latticci^ When MnAs is deposited on a GaAs substrate as 
a thin film, MnAs undergoes a strain induced increase in 
Tc and moreover develops a two phase region where the 
a- and /3-phases coexist between approximately 280 and 
340 K4 The magnetic transition temperature depends on 
choice of substrate, e.g. MnAs/GaAs(001) has a Tc=320 
K while MnAs/GaAs(lll)B has a Tc=346 as weU 
as on strain, e.g. stretching of the c-axis and compres- 
sion of the a-axis induce a reduction of Tci^ Even though 
the films arc fully relaxed after just a few monolayers^ 
the strain induced Tc increase seems to prevail in thicker 
films and through application of high magnetic fields the 
/3-phase can be transformed back to the a-phase even 
above Tc due to the structural strain.— The magnetic na- 
ture of the /3-phase has been under debate where some 
consider it to be paramagnetic^ while later suggestions 



point out that it has an antifcrromagnctic (AFM) short 
ranged structure?^ albeit without long range order In- 
vestigations of the magnetic anisotropy have determined, 
counter intuitively, that for MnAs the easy axis of mag- 
netization lies along the a-axis ([1120]-direction) and the 
hard axis along the c-axis ([0001]-direction).— Since both 
c and a lie in the film plane of MnAs/GaAs(001) films the 
magnetocrystalline anisotropy can be measured by sim- 
ply rotating the film 90°. An additional contribution to 
the in-plane anisotropy is shape-anisotropy arising from 
the striped structure with alternating a- and /3-stripcs, 
stretching out along [0001], that decrease and grow in 
width, respectively, as the temperature increases.— Re- 
ported values of the anisotropy constant^ are close to 
values obtained for bulk MnAs (7.4 x 10^ J/m"^ at room 
temperature) ^ 

We present an investigation of MnAs films grown on 
GaAs(lOO) substrates and describe how the magnetic 
properties depend on strain and uniaxial structure of the 
thin film. The temperature dependent 1st and 2nd order 
terms {Ki and K2) of the magnetocrystalline anisotropy 
(MCA) are derived from the free energy density equa- 
tion revealing an easy plane anisotropy in MnAs films 
at all temperatures. X-ray magnetic circular dichroism 
(XMCD) measurements reveal that the orbital moment 
does not strictly follow the spin magnetic moment in the 
temperature range where both the a- and /3-phases co- 
exist. 
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II. EXPERIMENTAL 

The MnAs films were grown in a KRYOVAK MBE 
system, using epi-ready GaAs wafers with (100) orien- 
tation as substrates. The substrates were attached to 
molybdenum holders by liquid In gluing, which ensures 
uniformity of the substrate temperature during the MBE 
growth. After introduction into the MBE growth cham- 
ber, the substrates were subjected to the typical proce- 
dures of thermal evaporation of native oxides and high 
temperature growth of a GaAs buffer layer. An As 
cracking cell (by DCA Instruments) operated at 900° C 
was used as the arsenic source. After the buffer layer 
growth the temperature of the As cracker was lowered 
to 600° C, which changes the As flux from As dimer rich 
to As tetramer (AS4) rich and the substrate tempera- 
ture was lowered to about 250 °C. This temperature was 
measured by an IR pyrometer operating in the 100 - 
700°C range. During the cooling time the substrates 
were exposed to the As4 flux. The desired surface re- 
construction (c(4x4) or d(4x4)) of the GaAs substrate 
was obtained by following established proceduresJ^ii^ 
The MnAs growth was started after reaching and sta- 
bilizing the substrate temperature at 250°C; the MnAs 
growth rate was chosen to be about 200 A/h. The growth 
was monitored by a reflection energy electron diffraction 
(RHEED) system. The RHEED images evolved rather 
quickly from two-dimensional (2D) 4x4 streaky patterns 
originating from the GaAs surface to diffraction spots 
indicating formation of MnAs islands during the initial 
growth stage. When the growth proceeded the RHEED 
images from MnAs:GaAs(100) evolved to an admixture of 
three-dimensional spotty and 2D streaky patterns, which 
were visible up to a MnAs film thickness of about 500 A. 
The spotty patterns eventually disappeared for thicker 
MnAs films. To prevent MnAs surface oxidation after 
the MBE growth the samples were capped with a 1000 A 
thick amorphous As layer deposited in the MBE chamber 
after switching off the substrate heater and cooling down 
the sample to the lowest temperature possible to attain 
in the MBE system. 

One of the investigated samples had a c(4x4) surface 
reconstruction of the GaAs(lOO) surface prior to MnAs 
growth, henceforth referred to as the c(4x4) sample, 
while the As- and de-cap samples had a d(4x4) surface 
reconstruction of the GaAs(lOO) surface prior to MnAs 
growth. The de-cap sample was subjected to an addi- 
tional post heat treatment at 310°C for 60 minutes (tak- 
ing place during removal of the As capping layer) in vac- 
uum. 

X-ray diffraction (XRD) 9 — 29 scans were carried out 
on a Bruker D8 Bragg-Brentano system, equipped with 
a CuXqj monochromator and a Vantec detector. Pole 
figures were measured using a Phillips MRD system with 
point focus and an X-ray poly capillary lens and a sec- 
ondary 0.18° parallel plate collimator and a graphite 
monochromator. This set-up was also used to measure 
the c-axis ((0004) peaks) of the MnAs (lOTO) orientation 



using grazing incidence XRD, with an incoming angle 
of about 0.7-0.8°. For the MnAs (lOTl) orientation an 
asymmetric scan was used to measure the c-axis. The 
MRD system was also used for rocking curve measure- 
ments, using a 2-bounce hybrid monoehromator/mirror 
and secondary X-ray mirror set-up. Due to the differ- 
ence in scattering power of the (1010), (1011) and (1012) 
planes the peak intensities in the 6-29 diffractogram 
shown in Fig. [T] are not proportional to the amount 
of each orientation. The (1011) orientation is overrepre- 
sented due the high scattering factor of this plane (about 
40 X higher than for the (1010) plane). The proportion 
of each orientation was therefore estimated by comparing 
the intensities in (1120) pole figures, cf. Fig. [2l 
A transmission electron microscopy (TEM) study was 
performed as a function of temperature between 89 and 
325K on a d(4x4) surface reconstructed sample using 
a JEOL-2100 microscope operated at 200keV with a 
double-tih cooling holder (Gatan CHDT3504). The 
cross-section TEM specimen was obtained, after mechan- 
ical lapping and polishing, by ion-milling (JEOL lon- 
slicer) thinning. The selected area electron diffraction 
(SAED) patterns along [101] GaAs (or [lOTO] MnAs), 
which were recorded on a CCD camera (Gatan ES500W), 
were used to determine the a- and c-axis at different tem- 
peratures. As a reference, the (0-20) GaAs refiection was 
also measured. Two selected areas of the sample, each 
with an effective size of 800 nm in diameter, were inves- 
tigated but no significant difference in the results were 
found between the studied areas. 

X-ray magnetic circular dichroism (XMCD) measure- 
ments were performed at beamline DlOll at the MAX- 
lab synchrotron facilityJ^ All measurements were per- 
formed by magnetizing the sample with an in-plane mag- 
netic field in the easy magnetization direction after which 
the absorption spectrum was measured in remanence us- 
ing total electron yield. The measurements were per- 
formed using 75±5% circularly polarized light and with 
the sample 13° from grazing incidence. The same pro- 
cedure was repeated after magnetizing the sample in the 
opposite direction. Saturation effects which are promi- 
nent at low angles are compensated for according to 
Nakajima et alJ^ The spectra are normalized at low pho- 
ton energies after which a step function is subtracted. 
We have used a correction factor for the spin magnetic 
moment of 1.47, which has been suggested to remedy 
problems arising when attempting to separate overlap- 
ping contributions (due to the relatively small spin-orbit 
splitting of the Mn L shell levels) from the L3 and L2 
edges 

Magnetization measurements were performed in a Quan- 
tum Design MPMS-XL Superconducting Quantum Inter- 
ference Device (SQUID) magnetometer. Magnetization 
(M) versus temperature (T) was studied between 5 K and 
400 K following a field cooled protocol. A weak magnetic 
field {H = 50 Oe) was applied along the crystallographic 
[1120], the easy axis of magnetization, at 400 K and the 
magnetization was measured as the sample cooled down 
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to 5 K. Magnetization versus field measurements were 
performed at different temperatures in the range 10 - 340 
K between 50 kOc and -50 kOe. The field was applied 
either along [1120] or along [0001], the in-plane easy and 
hard axis of magnetization, respectively. 
To corroborate the experimental studies we performed 
theoretical investigations based on density functional the- 
ory. In these studies, we have used linearized muffin tin 
orbital (LMTO) calculations with no shape approxima- 
tion for charge density or potentialJ^ The lattice is de- 
scribed by a non-overlapping spherical region surround- 
ing atomic sites called muffin-tin sphere and in between 
them an interstitial part. The basis set is formed with 
augmented linear muffin-tin orbitals. In-side the muffin- 
tin sphere basis function, charge density and potential are 
expanded in spherical harmonic functions up to Imax = 8 
along with radial functions and in a Fourier series in 
the interstitial region. In the interstitial region the basis 
functions are Bloch sums of spherical Hankel and Neu- 
mann functions. The Dirac equation is solved for a core 
charge density which means that no frozen core approxi- 
mation is considered. A multi basis formalism is used to 
ensure that all wave functions are well converged, i.e. we 
have considered three 4s, three 4p and two 3d orbitals 
for Mn and three 4s, three 4p and two 4d orbitals for As 
in the wave function. The modified tetrahedron method 
is used for integration in reciprocal space over 1152 k- 
points in the whole Brillouin zone (BZ) for self-consistent 
ground state calculations and to ensure convergence of 
the magnetic anisotropy energy (MAE), we carried out 
up to 9216 k-points in the full BZ. The exchange correla- 
tion potential is approximated by the PBE96 functional 
within the generalized gradient approximation. We also 
performed studies of the influence of electron correlations 
using dynamical mean field theory. 



III. RESULTS 
A. Structural properties 

Several out-of-plane and in-plane orientations have 
been reported for MnAs films deposited on GaAs (100) 
substrates. These orientations have been classified as A 
or B depending on the in-plane relationship with respect 
to the substratci^"— The A and B orientations have a 
twin relationship of 90° where in the ^-type the c-axis of 
MnAs is parallel with the GaAs [110] direction and for 
the S-type the MnAs c-axis is parallel to the GaAs [110] 
direction. Moreover, an index in the notation for the ori- 
entation was introduced to distinguish between in-plane 
and out-of-plane orientations of the MnAs c-axis . ^^'^^ For 
^0 and Bq there is an in-plane orientation of the MnAs 
c-axis, while for Ai, Bi and A2, B2 the c-axis is tilted 
out-of-plane with respect to the GaAs(OOl) surface (see 
Table U for details) . The orientation of the MnAs film 
can be controlled by the template layer. More in par- 
ticular, a very As-rich template with a d(4x4) surface 



reconstruction yields an Aq orientation while reduction 
of the As-concentration on the template yields a Bq or 
even a Bi orientation.— 
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FIG. 1: (Color online) X-ray diffraction 6-26 scans for the As- 
cap (red), de-cap (blue) and c(4 x 4) (black) samples obtained 
at 297 K. All samples show good epitaxy but with several 
orientations present in the films, see Table U] 



The d(4x4) films show a dominating (1010) orientation 
(Aq). Small amounts of additional (1011) and (1012) out- 
of-plane orientations were also observed. The fraction of 
(1011) orientation was estimated to about 5-13 %, while 
the (1012) orientation was less abundant, 0.4-2%. The 
de-cap film shows less (1011) and (1012) out-of-plane ori- 
entations, probably due to the heat treatment during the 
As-layer removal. The Aq orientation exhibits a rock- 
ing curve Full Width Half Maximum (FWHM) value of 
about 0.5° for the de-cap sample and 0.65° for the As- 
cap sample, while the Bi orientations have rocking curve 
FWHM values of 0.4-0.5°. The B2 orientation shows an 
even higher FWHM value, 1.2° and 1.9° for the As-cap 
and de-cap samples, respectively. These relatively high 
rocking curve FWHM values indicate a high degree of 
strain relaxation in the a- and 6-axis directions. 
The c(4x4) film exhibits a dominating (1011) orienta- 
tion with smaller amounts of the (1012) and (2023) ori- 
entations, 2% and 25%, respectively. The rocking curve 
FWHM value for the dominating orientation is about 
0.5°, while the (1012) orientation exhibits a higher value, 
about 1.0°. The (2023) orientation exhibits two rocking 
curve peaks, ±1.7° from the 6 angle. The FWHM value 
of the (2023) rocking curve peak is about 1.5-1.6°. This 
indicates that the fraction of (2023) orientation could be 
underestimated due to the high FWHM value. The films 
are almost completely relaxed with a exhibiting values 
close to the bulk value, while a strain up to 1% is ob- 
served for c (see Table |n| . In order to investigate the 
in-plane orientation and the fraction of each orientation 
present in the films, pole figures were measured. From 
the (1120) pole figure, shown in Fig. [2] for the de-cap 
and c(4x4) samples, the relative fraction (see Table HI)) 
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TABLE I: Epitaxial orientations of MnAs films de pos ited on 
GaAs (fOO), using the notation defined in Refs. [2ll - [24l . as 
well as the corresponding out-of-plane angles (a) of the c-axis 
with respect to the GaAs (001) surface. 



Orientation (hkl) 


Miller-Bravais 


a 


Notation 


(100) 


1010 


0° 


Ao, 


Bo 


(101) 


1011 


29° 


Ai, 




(102) 


10T2 


41° 


A2, 


B2 


(203) 


2023 


50° 


As, 


Bs 



of each orientation could be estimated through the height 
of each peak. It should be noted that A/S-orientations 
with the same a will appear on the same ■(/'-angle in the 
(1120) pole figure. 
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FIG. 3: Temperature dependence of the c/a-ratio. The c- and 
a-axis data were extracted from SAED patterns. 



(0.7%) of a. The results also indicate that the increase 
of c starts at slightly lower temperature compared to the 
decrease of a. 



FIG. 2: 1120 pole figures in linear scale with indications of 
orientations {Ao, B\, B2 and Ba) for a) the de-cap sample 
with mostly Ao orientation and a very small Bi contribution 
and for b) the c(4x4) sample with clear Bi and Bs orientation 
peaks. 

In the c(4x4) sample, the out-of-plane orientations [Bi 
and By,) dominate and almost no in- plane orientations 
of c {Ao/Bq) with respect to the substrate surface are 
found. For the d(4x4) surface reconstructed samples the 
Aq orientation dominates and through heat treatment a 
relaxation of the c-axis occurs (0.3%) together with an 
increase of the fraction of Ao orientation. Both the de- 
cap and As-cap samples have a fraction of i3-orientation, 
which is responsible for the small but finite value of the 
remanent magnetization when measuring the magnetiza- 
tion with the field applied along the hard magnetization 
direction, cf. Fig. [51 The out-of-plane A-orientations (Ai 
and A2) will also decrease the remanent magnetization 
when applying the field along the easy axis of magneti- 
zation (to be discussed below). 

Figure 13] shows the temperature dependence of the c/a- 
ratio as derived from the SAED patterns. Even though 
the SAED patterns are dominated by the a-phase reflec- 
tions, some weak reflections due to the /3-phase can be re- 
solved at the highest temperatures. The c/a-ratio starts 
to increase at about 250 K, which is where we expect the 
low temperature limit of the a — (3 coexistence region to 
be (see below), and reaches a value close to 1.55 at the 
highest temperature. This increase of the c/a-ratio is due 
to a 1% increase of c and a somewhat smaller decrease 



B. Magnetic properties 

From the temperature dependent magnetization curves 
shown in Fig. 2] it is seen that the highest Tc is obtained 
for the c(4x4) sample followed by the de-cap and As- 
cap samples (see Table HIl for values). For comparison, 
the temperature dependence of the remanent magnetiza- 
tion for the de-cap sample (M^) along with the magnetic 
moment obtained from XMCD measurements (rnxMCo) 
normalized with respect to Mr are plotted in Fig.U] The 
archetypical behavior of the field dependent magnetiza- 
tion in MnAs/GaAs(100) thin films at different temper- 
atures and fields applied along [1120] and [0001], here 
represented by the results obtained for the de-cap sam- 
ple, is seen in Fig. [S] When measured along the easy 
axis of magnetization (closed symbols) the film shows 
a very square like and narrow hysteresis curve, with a 
coercivity around 160 Oe at 10 K. The hard axis mag- 
netization exhibits a nearly linear relationship with the 
applied field up to saturation with a saturating field of 
« 35 kOe at 10 K. That the hard axis magnetization ex- 
hibits a small but finite remanent magnetization as well 
as magnetic hysteresis is related to the small fraction 
of out-of-plane B orientation in the film, as seen in the 
XRD-measurements, i.e. a small part of the film will 
have easy axis or close to easy axis alignment. 

MnAs has its easy axis of magnetization along the 
a-axis ([1120]-direction) and the hard axis along the c- 
axis ([0001]-direction) in the hexagonal unit cell. When 
MnAs is deposited on GaAs(lOO) both the a- and c-axis 
will be oriented in the fllm plane and thus the magnetic 
anisotropy can be studied by performing M versus H 
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FIG. 4: (Color online) Magnetization (M) versus tempera- 
ture (T) for the de-cap (blue squares), As-cap (red circles), 
and c(4 x 4) (black triangles) samples. The remanent magne- 
tization {Mr, purple diamonds) for the de-cap sample and 
the magnetic moment (mxMCD, green triangles) obtained 
from XMCD measurements are included for comparison. The 
TixMCD data have been normalized with respect to Mr at 100 
K. 

measurements along the in-plane directions defined by 
these two axes. In case of uniaxial crystal structure, the 
expression for the magnetic free energy density (F) is 
given by: 

F = Ki [flo cos^ uq + ai cos^ ai -t- . . . ] sin^ 9 + 
+ K2 [ao cos^ Q!o + fli cos^ «! + ...] sin"' 9 — 
- fioHMsCos9, (1) 

where Ki and K2 are the magnetic anisotropy constants, 
ao, ai. . . are the fractions oi Aq, Ai. . . orientations, re- 
spectively, in the sample and ao, ai. . . are the corre- 
sponding out-of-plane angles of the c-axis with respect to 
the film plane. The last term is the Zeeman contribution 
to the free energy and 9 is the angle between the magneti- 
zation and applied field directions. Through the relation- 
ships M = MgC0s9 and -^F = an expression for field 
dependence of the hard axis mag netization, M^'"''^/Af^ 
vs. H can be derived, which when fitted to the corre- 
sponding experimental data yields the anisotropy con- 
stants. The values for Ki and K2, given in Table HIl 
have been obtained by subtracting the easy-axis magne- 
tization contribution in the M^"'^'^ vs. H data before fit- 
ting the M^"'^'^ /Ms expression to the experimental data. 
As seen in Table HIl the As-cap and de-cap samples show 
a somewhat higher anisotropy and Tc than obtained for 
bulk MnAs,— while the c(4 x 4) sample displays lower 
iT-values but a higher Tc value. 

A difference in magnetic behavior can also be seen 
in the field dependence of the magnetization, see Fig. [S] 



FIG. 5: (Color online) Magnetization [M) versus field [H) at 
different temperatures in the range 10 — 305 K, for the de- 
cap sample with filled symbols corresponding to field applied 
along the easy axis ([1120]-direction) of magnetization and 
open symbols to field applied along the hard axis ([0001]- 
direction) of magnetization. 




H (kOe) 



FIG. 6; Magnetization [M) versus field {H) measured at 10 
K along the easy axis (filled symbols) and hard axis (open 
symbols) of magnetization, with the As-cap sample in red 
circles, de-cap sample in blue squares and c(4x4) sample in 
black triangles. 



where the closed and open symbols correspond to mea- 
surements along [1120] and [0001], respectively. A clear 
difference in the values for the saturation magnetization 
{Ms) is seen between the samples. We believe that the 
differences in Tc, Ms and MCA between samples are 
strain related and connected to the existence of multi- 
epitaxial orientations in the films. The c(4 x 4) sample 
exhibits a dominant out-of-planc orientation that is re- 
sponsible for the more rounded shape of the M'"""'* vs. 
H curve as well as for the higher Tc<^ Even though the 
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TABLE II: Anisotropy constants Ki and K2 measured at different temperatures (Tm), 7c, room temperature values for c- 
axis length and c/a-ratio, and fraction of each in- and out-of-plane orientation for the different samples of 1000 A thick 
MnAs/GaAs(100), compared with bulk values from Ref. IT^. 



sample 


T„, (K) 


_ft'i(xl05 J/m^) 


A"2(xl0-"' J/m^) 


Tc (K) 


Chcx (A) 


c/a 


Ao / Bi / A2 / B2 / B3 


As-cap 


10 


-14.0 


1.4 


317 


5.772 


1.555 


0.86 / 0.13 / 0.017 / - / - 




100 


-13.9 


1.3 












200 


-12.4 


1.2 












250 


-10.3 


0.9 












280 


-8.4 


0.9 












290 


-7.5 


0.9 












300 


-6.2 


0.6 












305 


-5.9 


0.4 










de-cap 


10 


-14.4 


1.7 


320 


5.755 


1.551 


0.95 / 0.050 / - / 0.004 / - 




100 


-14.4 


1.7 












200 


-12.8 


1.5 












250 


-10.6 


1.2 












280 


-8.8 


1.1 












290 


-7.9 


1.1 












300 


-6.5 


0.8 












305 


-5.8 


0.6 












1 n 


- 1 .0 


1.0 




^ 7Af{ 
0. ( 




U.UUo / U. / Z / - / U.UZ'l / u.zo 




100 


-5.6 


1.1 












200 


-4.8 


0.9 












300 


-3.2 


0.8 












320 


-2.9 


0.9 










Bulk 


77 


-12.0 




313 


5.712 


1.540 






198 


-11.0 














273 


-9.0 














299 


-7.6 













out-of-plane angles are compensated for when deriving 
Ki and K2, the obtained Ki values are a factor of 2 
smaller than the values obtained for the As-cap and de- 
cap samples as well as values derived for bulk samples 
It should be noted that the easy-axis magnetization mea- 
surements on the c(4x4) sample are made with the field 
applied along GaAs[110] instead of GaAs[110]. Although 
the As-cap and de-cap samples have a more uniaxial ori- 
entation with 86% and 92% of Ao-orientation, respec- 
tively, they still show some dissimilarities in Tc, as 
well as in Ki and K2 ■ This clearly corresponds well to the 
early investigations on bulk MnAs describing changes in 
magnetic properties with strain^ whereas in the present 
study the unit cell distortion comes from lattice mis- 
match and multi-orientation induced strain. The increase 
in Tc with the relaxing c is well in line with reports from 
likawa et al.^ who showed a relationship between com- 
pression/stretch of the lattice and changes in Tc. We 
would like to point out that the variations in Tc are due 
to two effects; the fraction and angle (a) of c-axis out-of- 
plane orientations^^ and the stretching/relaxation of the 
c-axis. The interdependence between a and Tc is evident 
for the c(4x4) sample where no in-plane alignment is de- 
tected. Although the c(4x4) sample has a close to bulk 
c/a-ratio (see Table HI)) the anisotropy seems to be more 
dependent on c-axis alignment along GaAs[110] than on 
c-axis strain. A conclusion that becomes even more ev- 
ident when examining the c/a-ratio for the As-cap and 



de-cap samples. 



1. X-ray Magnetic Circular Dichroism 

Mn ^2,3 x-ray magnetic circular dichroism spectra 
recorded at 100 K for the de-cap sample are shown in 
Fig. [7] together with the corresponding dichroic differ- 
ence. The spectra were obtained by reversing the mag- 
netization of the sample keeping the polarization fixed. 
Similar spectra were recorded for several temperatures in 
the interval 100-330 K. The spectra have been normal- 
ized to a common edge jump corresponding to a per atom 
normalization of the absorption cross section.— The use 
of the XMCD sum rulesS^ for Mn has been controversial 
since there is a mixing between the L2 and L3 edges and 
therefore a constant of 1.47 has been proposed to com- 
pensate for this effect Using a d-hole count of 4.9 as the 
number of empty d states in the sum rules, a magnetic 
moment of ^5 fi b /Mn atom is obtained at low temper- 
atures. This number is clearly too large as compared to 
the numbers obtained from magnetization measurements 
and theory. This discrepancy might have several expla- 
nations. We first note that the spin-quadrupole term (T^) 
was neglected in the sum rule analysis. A configura- 
tion would in a localized picture give the same occupation 
of all orbital levels resulting in a zero (T^) termi^ How- 
ever, one could of course argue that this approach might 
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not be strictly valid for the MnAs system due hybridiza- 
tion effeets and a deviation from a strict configuration. 
From our calculations we find that the d occupation is 
5.13 electrons. We have therefore also made explicit cal- 
culations of the {Tz) correction, following the recipe by 
Wu and Freeman.— In no case we find that the correc- 
tion is greater than 0.18 ^b, i-e., a correction less than 
5 percent of the total moment. We also note that with- 
out the corrective factor suggested in Ref. [H, we obtain 
a spin moment that is very close to the value obtained 
from calculations. Errors in the estimate of the number 
of holes, obtained from theory, can safely be ruled out, 
as this would require a charge transfer of almost 1.5 elec- 
trons not accounted for by the theoretical calculations. A 
more likely explanation for the discrepancy in magnetic 
moments connects to the value used for the correction 
factor and a more thorough investigation on the appli- 
cability of the corrective factor for the XMCD spin sum 
rule for MnAs is under preparationi^ To avoid distract- 
ing the focus from the more interesting point regarding 
the behavior of the spin and orbital moments as a func- 
tion of temperature, we have chosen to present moments 
obtained from XMCD sum rule analysis as magnetic mo- 
ments per hole. The results are summarized in Fig.|8l As 
shown in Fig. |5]3 the orbital moment (to/) behaves differ- 
ent from the spin magnetic moment (tos) given in Fig. [5^. 
As most evident from the TO//TOs-ratio, seen in Fig. [5J:, 
the orbital moment decreases more slowly than the spin 
magnetic moment in the temperature range where the a- 
and /3-phases coexist. This behavior is not strictly ex- 
cluded by any physical principle as we are dealing with 
a structurally driven first order magnetic transition. We 
note however that the result is an indication of structural 
changes in the a-phase in the coexistence region. This 
conjecture is confirmed by the temperature dependence 
of c/a shown in Fig.[3l the very similar temperature de- 
pendence of c/a and nii/ms is striking. Moreover, the 
observed variation of the orbital magnetic moment is in 
accordance with the theoretical calculations, discussed 
later in the text and seen in Fig. [9l where it is shown 
that an increased c/a-ratio can yield a larger orbital mo- 
ment anisotropy. 



C. Theoretical Calculations 

In this section, we will present ab initio calculations 
to address some behavior observed in the experiments 
discussed above. Let us start with the discussions on 
the magnetic anisotropy energy. This is a very delicate 
property to calculate, since a very dense sampling of the 
Brillouin-zone is needed in order to get accurate results. 
We have done careful convergence tests, with respect to 
the number of k-points, and our final calculation was 
done for 9216 k-points in the full Brillouin-zone. The 
final self-consistent result of the MAE is shown in Fig. [HI 
as a function of c/a-ratio. In Fig. [5] we also show the 
energy difference between a FM phase and an AFM one. 
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FIG. 7: Mn L2,3 XMCD spectra recorded at 100 K and their 
corresponding difference spectrum. Spectra were obtained by 
reversing the magnetization direction in the sample. 
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FIG. 8: a) Spin magnetic moment (ms), b) orbital magnetic 
moment {mi) and c) the ratio mi/ms- 
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FIG. 9: (Color online) Calculations for different c/a-ratio, 
where black symbols (red symbols) correspond to keeping the 
c-axis (a-axis) fixed while varying the a-axis (c-axis). a) dif- 
ference in total energies between ferromagnetic and antifer- 
romagnetic alignments of Mn spins in the [0001] -direction, 
b) magnetic anisotropy energy (MAE) as a function of c/a- 
ratio, c) orbital moment anisotropy (OMA) as a function of 
c/a-ratio on the left hand axis with filled symbols and solid 
line and spin magnetic moment as a function of c/a-ratio on 
the right hand axis with open symbols and dashed line. 






as a function of c/a-ratio. Note from the figure that 
we have made c/a changes while keeping either the c- 
axis or a-axis fixed at bulk experimental values (c=5.712 
A and a=3.71 A). The bulk experimental c/a-ratio is 
1.54, and for this geometry the FM phase is stable (sec 
Fig. in]) over the AFM phase with '-^O.l cV per formula 
unit. We also note from the theoretical results that the 
easy axis lies in the a6-plane, and that the c-axis is a hard 
magnetization direction. This result is in good agreement 
with experimental observations, and the size of the MAE 
(0.25 meV/f.u. «12xl0^ J/m'^) is rather substantial and 
well comparable to our experimental results, as can be 
seen by a comparison to data in Table as well as to 
bulk resultsji^ From Fig. [5^ it is seen that the AFM 



FIG. 10; (Color online) Spin-polarized band structures where 
black and red symbols indicate spin-up and spin-down bands, 
respectively, a) c/a=1.45 for fixed c and varying a, b) 
c/a=1.55 and c) c/a=1.45 for fixed a and varying c. 



phase is stable when the c/a-ratio is 1.65 for a fixed c- 
axis. A rather interesting feature is displayed in Fig. |9Jd 
in that the dependence on the calculated MAE is rather 
weak when keeping the a-axis fixed and varying the c- 
axis, whereas the dependence is very strong when the 
c-axis is fixed and the a-axis is varied. The microscopic 
reason for the MAE lies in the electronic structure^ and 
in order to analyze the behavior in Fig. [HI we show in 
Fig. [TU] the energy bands for different c/a-ratio. We also 



display energy bands when the c-axis is varied and the 
a-axis kept constant, and vice versa, for c/a-ratios equal 
to 1.45 and 1.55 (experimental value). It is well known 
that the contribution to the MAE from the spin-orbit 
coupling in perturbation theory is from occupied bands 
just below the Fermi level to unoccupied bands just above 
the Fermi level. Hence, it is the energy bands around the 
Fermi level one should inspect if one wants to investigate 
the MAE changes shown in Fig. |9l The energy bands 
in Fig. [TU] demonstrate an overall larger sensitivity when 
modifying the a-axis while keeping the c-axis fixed, which 
is consistent with the data in Fig. ^ This is particularly 
pronounced along the symmetry direction F — ?> A — >■ L — > 
M. Here entire segments of the spin-down Fermi surface 
are removed by the strain, since the spin-down band that 
for the unstrained case crosses the Fermi level, becomes 
completely unoccupied. This only happens for a strain 
where the c-axis is fixed. 
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FIG. 12: Spin-polarized DOS for (top) GGA with U=0 eV 
(top), (middle) DMFT with U=2 eV and (bottom) DMFT 
with U=4 eV calculations. 
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FIG. 11: (Color online) Magnetic anisotropy energy (MAE) 
vs. orbital moment anisotropy (OMA) for different c/a-ratio. 
The calculations were done by keeping either the a-axis (red 
circles) or the c-axis (black circles) fixed. 



The final theoretical result we would like to discuss is 
the relationship between MAE and the anisotropy of the 
orbital moment when measureed along an easy and hard 
magnetisation direction, the so called orbital moment 
anisotropy (OMA) shown in Fig. From perturbation 
theory it has been shown^ that there is a direct cou- 
pling between the MAE and the OMA. The derivation is 
made under the assumption that cross-terms in the spin- 
orbit interaction (i.e. ^(/+s_ -I- Ls+)) can be neglected 
compared to diagonal terms (i.e. ^IzSz)- For a material 
composed of one element this relationship holds well, and 
there is a direct proportionality between the MAE and 
OMA. However, a generalization of this theory to mate- 
rials composed of more than one element, shows a more 
complex behavior since both the MAE and OMA may 
have contributions from spin-orbit effects of all atomic 
species of the material. In short, this means that for a 



material like MnAs, the MAE of the magnetic moment, 
which is primarily located on the Mn atom, is also influ- 
enced by the spin-orbit interaction of the As atom. Since 
there are cross terms in atomic site index to both the 
MAE and OMA, which are large when there is strong hy- 
bridization of electron states centered on different atomic 
sites, there is no reason to expect a linear relationship 
between the MAE and OMA. Figure [TT1 shows a compar- 
ison between calculated values of the MAE and OMA, for 
different values of the c/a-ratio. Again, the calculations 
were done by keeping either the a-axis or the c-axis fixed. 
Note from the figure that the linear relationship between 
OMA and MAE is lost. Hence, for this material it is 
misleading to measure the OMA and from this value try 
to draw conclusions about the MAE. This has previously 
been shown experimentally for Au/Co thin filmsi^ 

The electronic structure of MnAs has been analyzed 
in the pasl,^"— using LSDA and/or GGA. A relevant 
question for this material concerns the influence of elec- 
tron correlations, which might be important for the elec- 
tronic structure and hence also for the magnetic proper- 
ties. Mn compounds often display electronic structures 
which have a correlated electronic d-shell, that may even 
become localized due to strong electron-electron interac- 
tion. Examples of this arc LSMO (i.e Lai-^rSr^MnOs)^ 
and MnOi^ Electron correlation has recently been dis- 
cussed even for one of the elemental forms of Mn, the 
7-phase4^ Hence, it is important to investigate this as- 
pect in MnAs, and for this reason we have compared the 
electronic structure and the spin and orbital moments 
of MnAs using GGA and dynamical mean field theory 
(DMFT - a recent method to take into account electron 
correlations) 4i In this calculation we have used the im- 
plementation by Grechnev et. al^ In Fig.[T2]we compare 
first the electronic structure from GGA calculations with 
DMFT calculations, using a Coulomb U (i.e. the strength 
of the electron-electron interaction) of for GGA and 2 
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and 4 eV for DMFT. As seen in Fig. [12] one can notice 
some changes in the electronic structure with the varia- 
tion of the value of U. The variation in the spin-up DOS 
is more prominent than the spin-down one. However, the 
resulting exchange splitting does not differ much as dis- 
cussed below. Next we turn our attention to the spin- 
(nis) and orbital (m;) moments. From the GGA calcula- 
tion we obtain for the Mn d-orbital a ms=3.12 /is/atom 
and a m; =0.026 /is/atom. From DMFT calculations 
for U=2 eV we obtain ms=3.12 //^/atom and m;=0.028 
/Lts/atom, while for U=4 eV we get ■ms=3.11 ^g/atom 
and TO/=0.017 ^s/atom. Our calculations thus show a 
rather weak influence on the spin and orbital moments 
as a function of increasing values of U, implying that 
electron correlations outside of what is included in GGA 
are not of major importance for MnAs. For the GGA 
calculation wc obtain a total moment of 3.13 /^s/atom, 
which should be compared to earlier theoretical results 
of 3.06 i_L B / atom^^^ and to the bulk experimental value 
of 3.45 ns/atomM 

IV. DISCUSSION 

The experimental Ms values for the As-cap and de-cap 
samples are at least 10-15% larger than what has been 
reported for thin filroii or bulk MnAs4^ Values similar 
to our experimental values were obtained by Sanvito 
and Hill^ from calculations of the effect of cell volume 
expansion where a fiB was obtained due to the 

reduction of p — d hybridization when then unit cell was 
stretched. The necessary increase of the cell volume 
required almost unphysical values of a and c, but the 
important observation is the large change in magnetic 
moment with lattice distortion. Our theoretical results 
also show that the magnetic properties depend strongly 
on lattice strain, especially on variations of a, see Fig. [3};. 
Since XRD is a global probe, the measured a and c 
correspond to averages over the sample volume and thus 
locally the tensile strain could be much larger. In fact, 
our SAED results, where the size of the selected area 
is 800 nm, indicate are more strained a-axis. However, 
we refrain from depending solely on the TEM results 
since thinning during sample preparation could induce 
additional lattice strain. Still, changes in Mf, and 
with annealing show that even small relaxations of the 



structure play a significant role. 

Our calculated magnetic anisotropy energy and easy axis 
direction are in good agreement with the experiments. 
Moreover, the theoretical results suggest that the change 
in magnetic properties is more sensitive by straining the 
in-plane axis than the out-of-plane one, an effect which 
has been attributed to the change in band structures 
due to strain. The calculated magnetic moments are 
observed to vary significantly upon straining the system. 
Following this observation, one may argue that the 
experimentally measured large moments may arise from 
the locally strained regions of the samples. In short, 
our theoretical and experimental results are in line with 
each other. 



V. CONCLUSIONS 

In conclusion, our findings indicate a more complex 
relationship between strain and magnetic properties in 
MnAs films than a simple stretch/compression of the unit 
cell axes. A strong dependence of Ki, K2 and Tc on the 
angle and fraction of out-of-plane orientation has been 
established. The ab initio calculations demonstrate a 
strong dependence of the magnetic properties on the vari- 
ation of the in-plane axis. We believe that local structural 
strains that distort the lattice give rise to the enhanced 
Ms and MCA measured in our samples. The possibility 
to alter the magnetic properties through controlled film 
growth with maintained epitaxy gives MnAs strong po- 
tential as a material for various spintronic applications. 
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